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Method of and Apparatus for Computation 
of Unbiased Power Delay Profile 

background of the invention 
cross-reference to related applications 

This patent application incorporates by reference the entire disclosure of a U.S. Patent 
Application entitled Method of and Apparatus for Path-searcher Window Positioning, bearing 
5 Attorney Docket No. 53807-00075USPT and filed on the same date as this patent application. 
BACKGROUND 
Technical Field 

The invention relates to digital wireless communications systems in which delays of 
individual multi-path components of a time-varying fading channel are to be estimated. The 

10 digital wireless communications systems may include, for example, systems using Code 
Division Multiple Access (CDMA) RAKE receivers. The invention relates in particular to 
improving robustness of detection of new multi-path components in a radio propagation 
channel, as well as of tracking known paths, by reliably indicating a proper position for a path 
searcher window. 

1 5 History of Related Art 

In wireless communications, a physical channel between a transmitter and a receiver 
is formed via a radio link. In most cases, no antenna of the transmitter is narrowly focused 
towards the receiver. In addition to a possible direct path, many other propagation paths 
often exist between the transmitter and the receiver. The other propagation paths typically 

20 result from reflections from objects near the transmitter or the receiver. Rays with similar 
propagation distances combine at the receiver, depending on an instantaneous phase 
relationship, and form a distinct multi-path component. The effect of a combination of the 
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rays depends on the instantaneous phase relationship of a carrier wavelength and also on 
distance differences among the rays. In the case of destructive interference, the combination 
of the rays leads to a significant decrease in path-gain magnitude (i.e., fading). 

Performance of a CDMA receiver is improved if signal energy carried by many multi- 
5 path components is utilized. A desired improvement in CDMA receiver performance may be 
achieved via a RAKE receiver. In the RAKE receiver, each of a plurality of multi-path 
components is assigned a despreader (i.e., RAKE finger). Each of the plurality of 
despreaders is assigned a reference copy of a spreading code. Each of the spreading-code 
reference copies is delayed in time by an amount equal to a path delay of a corresponding 

10 multi-path component. Outputs of the respective despreaders are then coherently combined 
via a RAKE combiner to produce a symbol estimate. 

The RAKE receiver preferably uses knowledge of the multi-path delays and channel- 
impulse values for all detected paths. To achieve a best possible signal-to-noise ratio at an 
output of the RAKE combiner, signal energy from as many physical paths as possible should 

15 be collected. In addition, tracking as many different physical paths as possible (i.e., maximal 
utilized diversity) significantly improves signal-reception robustness, since the probability of 
a simultaneous deep fade of all paths is reduced. Simultaneous deep fade of all paths is a 
phenomenon that typically leads to serious block-error-rate (BLER) degradation. 

A propagation channel structure (i.e., absolute and relative delays of the individual 

20 multi-path components) does not typically remain constant over time. Due to relative 
movement of the transmitter, the receiver, and nearby objects, delays of existing paths may 
change, old paths may disappear, and new paths may appear. In addition, a frequency offset 
between respective circuits of the transmitter and the receiver gives rise to a clock drift. The 
clock drift generally manifests itself as a gradual time-axis movement of the entire delay 

25 profile. To ensure proper operation of the RAKE receiver, the changing delays of all known 
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multi-path components should be tracked and new paths should be discovered quickly after 
the new paths appear. 

Due to the physical channel structure, in most cases relative positions of the nearby 
objects change. Thus, path lengths of the new paths usually do not differ significantly from 
5 path lengths of the existing paths. The macro-structure of the channel (e.g., mountains or 
groups of buildings that cause signal reflections) changes relatively rarely. Therefore, most 
often, the delays of the new paths are relatively similar to those of the existing, known, paths. 
Therefore, the delays of the new paths may be detected by searching in a delay domain near 
the known delays of the existing paths. 

10 FIGURE 1 is a block diagram of a typical RAKE receiver. A RAKE receiver 100 

includes a delay estimator block 102, a channel estimator block 104, and a RAKE 
despreader/combiner block 106. Received data are fed to the delay estimator block 102. The 
delay estimator block 102 evaluates an impulse response of a channel over a range of 
possible delays of the channel. A resulting delay profile, which may be a complex delay 

1 5 profile or a power delay profile, may then be subjected to peak detection and detected peak 
locations reported to the RAKE despreader/combiner block 106 as delay estimates for the 
multi-path components. The delay estimates are also used by the channel estimator block 
104 to estimate corresponding complex channel coefficients by despreading a pilot sequence 
and possibly filtering results over time to reduce the effects of noise and interference. 

20 Channel parameters are estimated in collaboration between the delay estimator block 102, 
which determines temporal alignment of a despreader portion of the RAKE 
despreader/combiner block 106, and the channel estimator block 104, which estimates the 
complex coefficients to be used by a combiner portion of the RAKE despreader/combiner 
block 106. A noise-plus-interference power estimate is also made. 
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A simple approach to delay estimation involves evaluating an impulse response of a 
channel over an entire range of possible delays (i.e., maximal assumed delay spread) of the 
channel. A resulting complex delay profile or power delay profile may then be subjected to 
peak detection and detected peak locations reported by the delay estimator block 102 to the 
5 channel estimator block 104 and the RAKE despreader/combiner block 106 as delay 
estimates. However, processing and power-consumption expenses of frequent execution of a 
full path-searching routine are usually prohibitive. Therefore, typical implementations use 
path searchers with observation windows shorter than the full search area (i.e., a maximal 
assumed delay spread). In addition, for any practical delay estimation, a path search is 

10 periodically undertaken to re-scan the delay range with the purpose of detecting new paths. 

A delay-estimation algorithm employed by the delay estimator block 102 extracts the 
path positions and finds the power delays with sufficient accuracy once the path positions 
have been discovered by the path searcher. A path-searcher window is positioned so that new 
paths are included within the path-searcher window. Since it is known with sufficient 

15 probability that the new paths will appear in the vicinity, in terms of the paths' respective 
delays, of the currently-known paths, the path-searcher window is usually placed so as to 
cover the currently-known paths. 

An estimate g(ri) of a current power delay profile for delays zy (i e [1, M]) typically 
includes a set of recently-detected or currently-tracked paths, in which case the delays v x are 

20 usually not contiguous, gfrj may also represent a contiguous region (t,- = r 0 + iAr) over 
which the path search is conducted. Other ways of representing the power delay profile are 
also possible. 

A suitable start position / for a path searcher window of length N w needs to be 
determined. A typically-used method for determining a suitable path-searcher window start 
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position for a next path-searcher activation is based on computing a center of gravity (i.e., 
mean excess delay) of the presently-known power-delay-profile estimate. A center-of- 
gravity position estimate C is computed as follows: 

Given C, the path-searcher window is placed so that most of the channel power is covered by 
the window. Because of space loss, a typical shape of the power delay profile exhibits 
exponential decay, such that the energy is concentrated towards the beginning of the region 
10 of interest. For reasonable coverage, the window can be placed, for example, 1/3 ahead of 

and 2/3 behind the value of C (i.e., / = C - ^ ). 

3 

In the case of a compact true power delay profile and a high receiver signal-to-noise 
ratio, C gives a consistent and reliable estimate of the true energy concentration in the 
channel. However, when the energy in the channel is distributed over a wide delay spread 
1 5 and when the signal-to-noise ratio of the power delay profile is poor, C is not so reliable. The 
noise-induced component of g(r,) causes a bias term that shifts the result of C towards an 
average non-power-weighted delay of all entries of the power delay profile. The size of the 
bias term depends on how far from each other the true center of gravity and the average delay 
are separated and also depends on the signal-to-noise ratio. In many practical cases, the bias 
20 term is large enough to shift the path-searcher window away from significant portions of the 
true power delay profile. 

To counteract noise-induced bias effects, g(r.) may be thresholded, which removes a 
portion of the noise-only samples and reduces bias. However, efficient noise removal 
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assumes the use of a rather high threshold, which may also remove channel components 
from, and thus distort, the power delay profile. 

The noise effect may also be reduced by noise subtraction, where average noise power 
cr* in the power delay profile is estimated. Instead of g(r.) , g(r.) - a] is used in the 
5 center-of-gravity computation. In a typical implementation, the center-of-gravity 

computation is based on a coarse power-delay-profile estimate over a range of N w delay 
values. The positions of the N p largest peaks of that power delay profile are taken as the 
delay values r r The noise floor cr^ 2 is estimated by averaging the N n = N w -N p smallest 
power-delay-profile values. However, this approach severely underestimates cr/ . Since the 
1 0 noise floor is not removed completely, significant residual bias effect remains. To improve 
the robustness and precision of the center-of-gravity computation and of the resulting path- 
searcher window placement, an approach is needed that more adequately removes the noise 
floor from center-of-gravity computations. 

SUMMARY OF THE INVENTION 

1 5 Embodiments of the present invention provide a method of and apparatus for 

computation of an unbiased power delay profile. A method of determining a noise-corrected 
power delay profile includes determining a power delay profile and calculating a noise- 
corrected power delay profile. The step of calculating the noise-corrected power delay 
profile includes using a biased noise-floor power estimate, the power delay profile, and a 

20 noise-scaling factor. 

An apparatus for determining a noise-corrected power delay profile includes a 
channel estimator, a despreader, and a delay estimator. The delay estimator is interoperably 
connected to the channel estimator and the despreader. The delay estimator is for 
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determining a power delay profile and calculating a noise-corrected power delay profile. The 
step of calculating the noise-corrected power delay profile includes using a biased noise-floor 
power estimate, the power delay profile, and a noise-scaling factor. 

An article of manufacture for determining a noise-corrected power delay profile 
5 includes at least one computer readable medium and processor instructions contained on the 
at least one computer readable medium. The processor instructions are configured to be 
readable from the at least one computer readable medium by at least one processor and 
thereby cause the at least one processor to operate as to determine a power delay profile and 
calculate a noise-corrected power delay profile. The calculation of the noise-corrected power 
10 delay profile includes using a biased noise-floor power estimate, the power delay profile, and 
a noise-scaling factor. 

BRIEF DESCRIPTION OF THE DRAWINGS 

A more complete understanding of exemplary embodiments of the present invention 
can be achieved by reference to the following Detailed Description of Exemplary 
1 5 Embodiments of the Invention when taken in conjunction with the accompanying Drawings, 
wherein: 

FIG. 1 , previously described, is a block diagram of a typical RAKE receiver; 

FIG. 2 is a flow diagram that illustrates calculation of a center-of-gravity position 
estimate in accordance with principles of the present invention; and 
20 FIG. 3 is a graph that illustrates exemplary values of a noise-scaling factor in 

accordance with principles of the invention. 
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DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS OF THE 
INVENTION 

5 Embodiments of the invention permit an unbiased center-of-gravity position estimate 

to be calculated given a set of N w power-delay-profile values. The center-of-gravity position 
estimate is computed from a subset of the power-delay-profile values. The subset may be, for 
example, the N p largest peaks. A noise floor component is removed from the N p peaks before 
the center-of-gravity position estimate is computed. A noise floor power is computed by 

10 averaging a subset of the power-delay-profile values. The subset of the power-delay-profile 
values may be, for example, the N n smallest values. A noise-scaling factor y is introduced 
into the center-of-gravity position estimate. The noise-scaling factor, which converts the 
biased noise floor power into an unbiased value, is computed based on a probability 
distribution of the noise component and on N W9 N p , and N„. As a result, the computed center- 

1 5 of-gravity position estimate is substantially unbiased and the path-searcher window is 
positioned with the best likelihood of covering all the paths included in the power delay 
profile. FIG. 2 is a flow diagram that illustrates calculation of a center-of-gravity position 
estimate C in accordance with principles of the present invention. A flow 200 begins at step 
202, at which step a noise-scaling factor y is determined. The noise-scaling factor y may 

20 be determined offline and utilized online by the delay estimator block 102 in determining the 
center-of-gravity position estimate C. At step 204, a power delay profile g m is determined. 
At step 206, a noise -corrected power delay profile is calculated using a biased noise-floor 
power estimate a g 2 , the noise-sealing factor y , and the power delay profile g m . At step 208, 

the center-of-gravity position estimate C is calculated using the noise-corrected power delay 
25 profile. 
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A path searcher processes a window of length N w (e.g., chips) and, at each nth 
activation, the power delay profile g m is estimated, with m e [m 0 ,m 0 + N w -l]. Out of the 
power delay profile, N p peaks (specifically delays r k and powers /** for k = 1 ... N p ) are 
detected and reported. Let A* = pk + where pk is the actual path power after power-delay- 
5 profile accumulation and tik is a residual noise plus an interference component. The power- 
delay-profile noise-floor power estimate based on the N n lowest power-delay-profile 
samples is also available. 

Traditionally, the center-of-gravity position estimate C is computed based on the 
reported delays and noise-corrected powers as indicated in equation (2) below. 

The noise floor power Z = cr g 2 is estimated as the mean of the smallest power-delay- 
profile samples. However, the noise included in a selected peak set {h k }, k = 1. . . N p , does 
not have the same distribution as the smallest N n samples. In the most troublesome cases, 
when there are many noise samples and few or no true paths, {hk} contains all the largest 
1 5 noise samples whose mean power (denoted by Z*) is significantly larger than the average 
power of the remaining noise floor. The are well-modeled as magnitude squares of 

complex Gaussian random variables. In other words, the /** have the X ^ -distribution with the 

.X 

1 

probability density function f(x) = . e g . We can thus express 

r xf(x)dx 
Z = Jo 



20 (corresponding to the near- zero portion of the probability density function) and 
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Z N p /N w 

(corresponding to the tail of the probability density function), where the integration limits a 
and p are defined so that J°° f(x)dx = ^- and j° f(x)dx = . 

The noise-scaling factor y is then determined as the following ratio: 

5 P=-§-' (3) 

AT 

FIG. 3 is a graph that illustrates exemplary values of y as a function of — - , assuming a 
N 

typical value of — - of 0.8. It is apparent from FIG. 3 that, in the example shown therein, 

using Z directly for noise correction underestimates the relevant noise power Z* by a factor 
of approximately 4-8. In other words, only approximately 12-25% of the noise power is 
1 0 removed. 

It may be shown that the expected bias of the center-of-gravity value equals 



where C true is the actual center of gravity. Noise removal may be achieved by minimizing 
E[nk -yZ]. Bias would in fact be removed if the true yZ = En k could be used, Equation (3) 
1 5 representing a practical estimate of such a y . As shown, the bias is removed by 

embodiments of the invention when a noise-scaling factor y is introduced into the center-of- 
gravity calculation as follows: 



11 

DALLAS2 1014951vl 53807-00076USPT 



Patent Application 
53807-00076USPT 



P T k (h k -yZ) 



The noise-scaling factor y may be determined in advance, thus avoiding on-line 
computations. The denominator of equation (4) directly yields the unbiased signal power 
estimate P h wherein 

Pi=Zh-M P rZ (5) 

k=\ 

and equation (5) may be utilized in other parts of the path-searcher algorithm. 

N p 

The value of y based on a fixed argument — - may differ from the optimal value 

when there are many actual paths and fewer noise-only entries within the observation 

window. The optimal noise-scaling factor is 

j__ N paths + y{N p — N paths ) 
N p 

where N pat h S is the number of true paths among the N p selected peaks. When the value of y 
differs from the optimal value due to many actual paths and fewer noise-only entries, an over- 
correction results in residual bias in the opposite direction to remain. However, the ratio 



, and correspondingly any center-of-gravity bias, is kept small. In high-noise 



cases, practical tests indicate a good model match and removal of almost all of the bias. 

Embodiments of the invention permit the actual noise floor of the observed power 
delay profile included in the center-of-gravity computation to be determined, based on a noise 
floor estimate derived from the partial noise distribution. As a result, the calculated center- 
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of-gravity value is unbiased or has a significantly reduced bias, improving the robustness of 
the window placement. In addition, the instantaneous center-of-gravity value calculated in 
accordance with principles of the invention may be time averaged or combined with other 
center-of-gravity values from other time or space points without incurring any bias-related 
5 errors and without requiring heuristic selection steps. 

Embodiments of the invention are computationally efficient, in that typically only one 
table lookup and one multiplication, in addition to the conventional center-of-gravity 
computation, are required. Given the relatively-infrequent computation of the path-searcher 
window position, the computational premium is usually negligible. 

10 As will be appreciated by those having skill in the art, the invention may take many 

embodiments. For example, embodiments of the invention may be used when other 
combinations of the number of peaks and the number of noise samples are used (e.g., so that 
N p + N n p±N w .). In such cases, computation of the noise-scaling factor will then be 
numerically different; however, implementation details will be obvious to a person skilled in 

1 5 the art. In addition, principles of the invention may be employed in any applications in which 
the center-of-gravity needs to be computed. A complex delay profile may be used instead of 
a power delay profile, in which case the power of each delay profile element is found by 
multiplying the complex coefficient by its complex conjugate. Either an averaged or an 
instantaneous power delay profile may be used. 

20 Although embodiment(s) of the present invention have been illustrated in the 

accompanying Drawings and described in the foregoing Detailed Description, it will be 
understood that the present invention is not limited to the embodiment(s) disclosed, but is 
capable of numerous rearrangements, modifications, and substitutions without departing from 
the invention defined by the following claims. 
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